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We present explicit analytical solution for the problem of minimization of the function F(x, y) = 
^ ^ 3 j=i m iV ( x — x i) 2 (y ~ Vi) 2 ' we nnc ^ the coordinates of stationary point and the correspond- 
^ ing critical value as functions of {irij,Xj,yj}^ =1 . In addition, we also discuss inverse problem of 
<^ finding such values for m 1 ,m 2 ,m 3 for which the corresponding function F possesses a prescribed 
'. position of stationary point. 



o 
u 



1 Introduction 



Given the coordinates of three noncollinear points Pi = (xx,yi),P 2 — (^2, Z/2), P3 = (x 3 ,y 3 ) in the 
CJ. plane, find the coordinates of the point P* = (x*,y*) which gives a solution for the optimization 
problem 



mmF(x,y) for F(x,y) = > m j\ / ( x ~ x j) 2 + (y ~ Vj) 2 ■ (1-1) 
(N- 3=1 

^ Here m 1 ,m 2 ,m 3 are assumed to be real positive numbers and will be subsequently referred to as 
^ weights. 

O The stated problem in its particular case of equal weights m\ = m 2 = m 3 = 1 is known since 1643 
^ as (classical) Fermat-Torricelli problem. It has a unique solution which coincides either with one of 
the point Pi,P 2 ,P3 or with the so-called Fermat or Fermat-Torricelli point [21 H] of the triangle 

Generalization of the problem to the case of unequal weights was investigated since XIX century. 
03 \ This generalization is known under different names: Steiner problem, Weber problem, problem of 
railway junctional [31 [S], the three factory problem pj. The two last names were inspired by the 
optimal transportation problem like a following one. Let the cities Pi,P 2 and P3 be the sources of 
iron ore, coal and water respectively. In order to produce one tonne of steel one the steel works need 
mi tonnes of iron, m 2 tonnes of coal and m 3 tonnes of water. Assuming that the freight charge for 
tonne-kilometer is independent of the nature of the cargo, find the optimal position for steel works 
connected with Pi, P 2 , P3 via straight roads so as to minimize the transport costs. 

In the rest of the paper this problem will be referred to as the generalized Fermat-Torricelli 
problem. Existence and uniqueness of its solution is guaranteed by the following result [1]: 
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Theorem 1.1 Denote by a\, «2,«3 the corner angles of the triangle PiP 2 P^. If the conditions 



m\ < m\ + m\ + 2 m 2 m 3 cos «i, 
wij < + rrig + 2 m!m 3 cos «2, 
mg < m 2 + m 2 , + 2 m\m<i cos «3 



(1.2) 



are fulfilled then there exists a unique solution P* = (x*,y*) G R 2 for the generalized Fermat-Torricelli 
problem lying inside the triangle P 1 P 2 Ps. This point is a stationary point for the function F(x,y), 
i. e. a real solution of the system 



If any of the conditions is violated then F(x,y) attains its minimum value at the corresponding 
vertex of the triangle. 

Let us overview some approaches for finding the point P*. The first one is geometrical: the point 
is found as intersection point of special construction of lines or circles. For the equal weighted case, 
Torricelli proved that the circles circumscribing the equilateral triangles constructed on the sides of 
and outside the triangle P1P2P3 intersect in the point P*; for an alternative Simpson construction of 
P* see [5] . For the general (unequal weighted) case see [3j [8] . 

The second approach is based on mechanical modeo a horizontal board is drilled with the holes 
at the points Pi,P 2 ,P3 (or at the vertices of a triangle similar to P1P2P3) ; three strings are tied 
together in a knot at one end, the loose ends are passed through the holes and are attached to 
physical weights proportional to mi,m2,m^ respectively below the board. The equilibrium position 
of the knot yields the solution [3]. 

The third approach, based on gradient descent method, was originated in the paper [H]; further 
developments and comments can be found in [T| [TU]. 

The present paper is devoted to the fourth approach — analytical one. We look for explicit 
expressions for the coordinates of stationary point P* as functions of {rrij,Xj,yj}^ =1 . Although the 
existence of such a solution by radicals, i.e. in a finite number of operations like standard arithmetic 
ones and extraction of (positive integer) roots, is not questioned in any review article on the problem, 
we failed to find in publications the constructive and universal version of an algorithm even for the 
classical (i.e. equal weighted) case. 

3 Sometimes incorrectly called Polya's mechanical model. 




(1.3) 
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2 Algebra 

Theorem 2.1 Under the conditions U.S\) . the coordinates of stationary point (x*,y*) of the function 
F(x,y) are as follows: 

_ K1K2K3 ( xi x 2 x 3 \ _ K X K 2 K 3 ( yi V2_ , V3_\ (0 n 

* " ASad Ui K2 kJ ' V * ~ ASad Kk, + K 2 kJ [ ] 



•with 



Here 



F(a;*,x*) = min F(x,y) = Vd . 

d= ^m\K x + m\K 2 + m 2 3 K 3 ) (2.2) 
la 

2S(y + \ hi( r i2 + r n - 40 + ^2(^3 + rf 2 - r 2 13 ) + m 2 3 (rj 3 + r 2 23 - r 2 12 )] . (2.3) 



and 



rjt = \l fa ~ + (%• - Vi) 2 = \PjPt\ for {j,£} C {1,2,3} 



S = \x x y 2 + x 2 y 3 + x 3 y t - Xiy 3 - x 3 y 2 - x 2 y x \ ; (2.4) 



1 

a — - 
2 



and 



\J-m\ 


— m\ — m 3 + 2 m\m\ 


+ 2 mfm§ + 2 rri^m 2 ; 


(2.5) 


K x = 


(r 2 2 + rj 3 - r 2 3 )a + ( 


'm 2 , + m| — ml)S, 




K 2 = 


(r 23 + r h - r 2 13 )cr + ( 


[m\ + m| — m 2 ,)^, 


(2.6) 


K 3 = 


(r 2 3 + r h ~ r\ 2 )o- + ( 


'm 2 + m 2 — m 3 )S. 





Proof. Firstly, we established via direct computations the validity of the following equalities: 

K X K 2 + K X K 3 + K 2 K 3 = AaSd , (2.7) 

and a dual one for (12.21) : 

r\ 3 K x + r\ 3 K 2 + r 2 2J K" 3 = 2 Sd . (2.8) 
Secondly, let us deduce the following relationships: 



y/(x. - Xj y + - % -) 2 = ^ /or J 6 {1,2,3}. (2.9) 



2(7^ 

To prove (12. 9p for j = 1 we first represent X* and ?/* given by (12. ip with the aid of (I2.7P : 



Thus, 



AaSd 



AaSd 

2 



%2 X 3 

K 2 K 3 



+ 



K 2 K 3 



m_ _ m_ 

K 2 K 3 



(x 2 - xi) 2 + (y 2 - yi) 2 , (x 3 - xi) 2 + (y 3 - 



jc2 



+2 



(x 2 - Xi)(x 3 - Xi) + (y 2 ~ V\)(Vz ~ Vi) 
K 2 K 3 



K X K 2 K 3 
AaSd 

-2 



(J^p + r\ 3 K\ + (r 2 2 + r\ 3 - r 2 3 )ir 2 K 3 ] 



(4aSd) 2 



(AaSd) 2 



[(r 2 12 K 3 + r\ 3 K 2 )(K 2 + K 3 ) - r\ 3 K 2 K 3 \ 



[(2 Sd - rl 3 K 1 )(K 2 + K 3 ) - r\ 3 K 2 K 3 \ 



(AaSdy 



[2 Sd(K 2 + K 3 ) - r\ 3 (K x K 2 + K X K 3 + K 2 K C 



123 



K 2 9 ^dK 2 

' 1 [2 Sd(K 2 + K 3 )-A r 2 23 aSd] = [K 2 + K 3 -2 r 2 23 a] 



(AaSd) 



J2U K 2 



[2 mlS] 



(AaSd) 2 
m\K 2 



8Sda 2 l-"m~J 4 a 2 d ■ 
Similar arguments hold for j G {2,3} in (12. 9p . To complete the proof of those equalities it should 
be additionally verified that the values Ki, K 2 , K 3 are nonnegative. This will be done in the next 
section. 

Now let us prove the first statement of the theorem. Substitute (12.1 p into the left-hand side of 
the first equation of (11. 3p . The resulting expression can be represented with the aid of formulae (12. 7p 
and flUD as 



l l l 

K~ 1 + K~ 2 + K~ 3 



Xi 



x 2 



Ki + K 2 + K 3 



x 3 \ &m 



0. 



Ki K 2 K 3 

Similar arguments are valid for the second equation from (jl.3p . Finally compute F(x*,y*): 



3=1 



^ 2aVd 2a\[d 



□ 
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Tests 1 



Pi 

mi 

2 



P2 
m 2 

(1,1) 
3 



P 3 
m 3 

(5,1) 
4 



Vd 

(2^(4103 + 1833^) , ^(29523 -4481^)) 
« (3.9086,1.41 52) 

W = 2v / 79+l57l5 « 23.4174 



/ 751 647 \ 

V 485 ' 485 / 

« (1.5484, 1.3340) 
-s/d = -v/970 « 31.144J 



2. 



(2,6) 
3 



(1,1) 
5 



(5,1) 
4 



3. 



(0,0) 
3/2 



(2,0) 
2 



■y/2, y/2) 
2 



1 



v / TlO ' %/2 v 7 ^ v 7 !!^ 

(0.0068,0.0165) 



Vd 



'32 + ^ + 3Jf « 7.9997 



(39,57) 
18 



(22,42) 
41 



(42,75) 
52 



« (37.0432,61.4053) 
= a/3068047 + SgiSv/TSlT w 1845.8994 



The exact coordinates of P* in Test 1.4 are as follows: 



296577529815837 , 357441196078431 /7F1T 
9297789607234 ' 6020318770684015 V 



271001243105952 , 432306390086253 /7F1 1 
4648894803617 ~ l ~ 12040637541368030 V ' 



3 Geometry 

Let us give an interpretation for some constants appeared in Theorem 12.11 Being rewritten in an 
alternative form, the constant (E 



S 



1 1 1 

det I x\ x 2 x 3 

yi V2 vz 



is recognized as the doubled area of the triangle PiP 2 P 3 . As for the constant ( 12. 5p . factorization of 
the radicand in its right-hand side leads one to the form 



a 



' mi + ^2 + m>3 ( ™>i + ^2 + 1^3 



mi 



mi + m 2 + m 3 



m 2 



mi + m 2 + m 3 



m 3 
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Figure 1: 

which can be treated as the Heron's formula for the doubled area of a triangle formed by the triple of 
weights mi,m 2 ,m 3 . Under the restrictions (11. 2ft . such a triangle exists. Construct now this triangle 
and denote its angles as shown in Figure 1. 

Then the first formula from (12.61) can be represented with the aid of the law of cosines as 

/ r{ 2 + rf 3 - r 23 m\ + m 2 3 - m\\ (2 r 12 r 13 cos a x 2 m 2 m 3 cos 0i \ 
Kl = aS { S + ; ) =aS { S + a ) 

= 2 aS(cot aix + c °t Pi) ■ 

On rewriting the first condition from ( II. 2p in the form cos«i + cos/?i > 0, one can conclude 
that cot oil + cot 0i > and, thus, K x > 0. In a similar way the expressions for K 2 and K 3 can be 
deduced, and established that under the restrictions (11.2j) they both are positive. This completes the 
proof of Theorem 12.11 

Remark. Set the dual generalized Fermat-Torricelli problem: let the triangle be composed of 
the sides with the lengths equal to mi, m 2 and m 3 ; let the weights r 12 , r 23 , r 13 be placed in its vertices 
as shown in Figure 2. 

The minimum value for the objective function will be the same as in the direct problem since (12. 3p 
is equivalent to 

oCi^r2/2, 2 2\,2/2, 2 2\ , 2 / 2 , 2 2\1 

2 Sa + - [r 12 {m 1 +m 2 - m 3 ) + r^m^ +m 3 - m 2 ) + r 23 (m 2 + m 3 - mjj . 
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Figure 2: 

4 Classical Fermat-Torricelli problem 

Consider now the case m\ = m 2 = m.3 = 1. 

Theorem 4.1 Let all the angles of the triangle P1P2P3 be less than 2n/3, or, equivalently: 

r i2 + r i3 + r i2^i3 - rj 3 > 0, r-23 + r\ 2 + ri 2 r 23 - r 2 3 > 0, r 2 3 + r 23 + ri 3 r 23 - r 2 2 > . 
T/ie coordinates of Fermat-Torricelli point for this triangle are as follows: 

X *~2V3Sd\h k 2 k 3 J ' V * ~ 2V3Sd \h k 2 k 3 J 

with the corresponding minimum value of the objective function: 

3 



Here 



and 



F(x*, £*) = min J (x — Xj) 2 + (y — y 3 ) 2 = \fd . 



1 T* 2 I yi 2 I yj 2 

d = —i={ki + k 2 + fes) = 12 + I 3 + 23 + v^S (4.2) 
V 3 2 



kl = ^( r L + r ?3 - r 2 3 ) + 5 , ^2 = ( r 23 + r L - r l 2 3 ) + S , k 3 = -^( r 13 + r 23 ~ r n) + S. 

with the rest of the parameters coinciding with those from Theorem \2.1\ 



It turns out that expressions (14. ip . being represented as rational fractions with respect to {xj, Uj}^ =1 , 
can be reduced further to the form where denominators become "area free" [T2] : 

Corollary. Under the conditions of Theorem 14.11 the coordinates of Fermat-Torricelli point are 
as follows: 



Here 



x* = ~ j= Kx 1 + x 2 + x 3 )\S\ + \/3 (xir^ + x 2 rl 3 + x 3 r\ 2 ) 
2y 3 a 

1 1 1 

+3sgn(S') yi y 2 y 3 

x 2 x 3 + 1/21/3 xi%3 + ym %iX2 + ym 

V* = ^j^[{yi + V2 + V3)\S\ + V3 (yirl 3 + y 2 rl 3 + y 3 rl 2 ) 

1 1 1 

-3sgn(S') xi x 2 x 3 

x 2 x 3 + y 2 y 3 X1X3 + y ± y 3 x x x 2 + y x y 2 



(4.3) 



(4.4) 



1 1 1 

S = det ( xi x 2 x 3 

yi y2 2/3 

Remark. Result of the last corollary can be extended to the generalized Fermat-Torricelli prob- 
lem: numerators and denominators of the formulae (12. ip can be reduced by the common factor S. 
We do not present here the resulting expressions since they are inelegantly ponderous [12J. 



Tests 2 





Pi 


P2 


P3 




y* 


v 7 ^ 


1. 


(1,1) 


(3,5) 


(7,2) 


efy (1029 + 79^) 
« 3.3939 


^ (1053 + 647^) 
« 3.1639 


^41 + 22^ 
« 8.8941 


2. 


(1,2) 


(3,3) 


(4,1) 


1(15 + V3) 
« 2.7886 


|(3 + V3) 
« 2.3660 


^10 + 5^ 
« 4.3197 


3. 


(0,0) 


(399,0) 


/5005 9555^ 
v 38 ' 38 / 


21255 
133 

« 159.8120 


8580^ 
133 

» 111.7368 


784 


4. 


(0,0) 


(2,0) 


(0,1) 


13 + 39 V 

« 0.2545 


13 39 v 

« 0.3045 


^5 + 2^ 
« 2.9093 



Test 2.3 is generated from [TT], test 2.4 is taken from [9]. 
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5 Inverse problem 



Given the coordinates of the point P* = lying inside the triangle P1P2P3, find the values for 

the weights mi, 777.2,7713 with the aim for the corresponding function F(x,y) to posses a minimum 
point precisely at P*. 

Theorem 5.1 Let the vertices of the triangle P1P2P3 be counted counterclockwise. Then for the 
choice 



|P*Pi| 
the function 



1 1 1 

X* %2 X3 

y* V2 2/3 



m 



|P*P 2 | • 



1 1 1 

X\ X* X3 

yi y* 2/3 



|p*p 3 



111 

X\ X2 x* 

y\ y2 y* 



(5.1; 



3=1 

has its stationary point at P*. Provided that the latter is chosen inside the triangle P1P2P3 the values 
( 15. 1\) are all positive and 



P(x*,y*) = mm F{x,y) 



1 


1 


1 


1 










y* 


yi 


2/2 


2/3 



xl + yl xj + yf x 2 2 + y 2 2 x\ + yl 



(5.2) 



Proof. Substitute x = x* , y = y^ and the values ( 15 - 1 p into the left-hand side of the first equation 
from fTO|) : 



(5.3) 



Represent this combination of the third order determinants in the form of the fourth order determi- 
nant, namely 





1 


1 


1 




1 


1 


1 




1 


1 


1 










+ (x* - x 2 ) 








+ (x* - x 3 ) 


Xi 


%2 






y* 


2/2 


2/3 




yi 


y* 


2/3 




2/1 


2/2 


y* 



1 


1 


1 


1 




Xi 


^2 


%3 


2/* 


2/1 


2/2 


2/3 



(expansion by its last row coincides with (I5.3P ). Now add the second row to the last one: 

1111 

X* Xi x 2 x 3 

y* 2/1 2/2 2/3 

ry> /y rf /y 
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In this determinant the first row is proportional to the last one; therefore the determinant equals 
just zero. The second equality from (11.31) can be verified in a similar manner. 
Let us evaluate F(x*,y*): 









1 


1 


1 




- xi) 2 + (y* 


- yi) 2 } 






^3 








y* 


2/2 


2/3 








1 


1 


1 


+ [(£* 


- x 2 f + (y* 


- yi) 2 ] 


X\ 




^3 








yi 


y* 


2/3 








1 


1 


1 


+ [0* 


- X3) 2 + (y* 


- ysf] 














yi 


2/2 


2/* 



To prove the equality (15. 2p let us split it into the x-part and the ?/-part. First keep the x-terms in 
brackets of the previous formula: 



t x* 



1 1 1 

X* X2 X3 

y* 2/2 2/3 



+ (x* 



%2 



1 1 1 

%\ 3^3 

2/1 y* 2/3 



+ (x* - x 3 y 



1 1 1 

X\ X2 x* 

2/1 2/2 y* 



Similar to the proof of the first part of the theorem, represent this linear combination as the deter- 
minant of the fourth order 



1 


1 


1 


1 






x 2 


£3 


y* 


2/1 


2/2 


2/3 



(x* - xi) 2 (2* - x 2 ) 2 (x* - x 3 ) 5 
Multiply the first row by (— x^), the second one by 2 a;* and add the obtained rows to the last one: 



1111 

X* X\ X2 £3 

y* 2/1 2/2 2/3 

2222 



(5.4) 



The y-part of the equality ( 15. 2 p can be proved in exactly the same manner with the resulting deter- 
minant differing from (15 .4p only in its last row. The linear property of the determinant with respect 
to its rows completes the proof of (15. 2p . □ 
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Remark. Solution of the inverse problem is determined up to a common positive multiplier, i.e. 
the solution triple (mi, m 2 , m 3 ) is defined by the value of the ratio mi : m 2 : m 3 . Up to this remark@, 
solution of the inverse problem is unique: we have proved this statement via direct computations 
starting from the formulae (12. ip . 

Example 5.1 Let P x = (2, 6), P 2 = (1, 1), P 3 = (5, 1) and 

P* = ( —J— (4103 + 1833^15^ , — — ( 29523 - 4481^ 
V2866 V / ' 8598 V 

Find the values for the weights m^m^m^ from Theorem \5.1\ 
Solution. Formulae (15. ip give: 

2(20925 - 4481^) 



18481401 




2(15105 - 2342 v 


/15 V' 


6160467 




8(— 1185 + 1598^ 




18481401 



316380606 + 35999826V15 



m* 2 = — — — - — ^1/75400161 - 9169767715, 



mt = — — -V 8335761 - 2050623^, 

3 18481401 v 

with 

F{x*,y*) = ——(-333980 + 193436^) . 
4299 

Now compare the obtained result with the one represented in Test 1.1: according with the last remark 
one might expect that 

m* : m* 2 : m* z = 2:3:4. 

Let us verify this fact: 



m* 1 20925 -4481 V15 316380606 + 35999826V15 



m* 2 3 15105 - 2342V15 V 75400161 - 9169767V15 

1 (20925 - 4481 y/TE) (15105 + 2342^) 
"3 145886565 



X' 



(316380606 + 35999826V15) (75400161 + 9169767V15) 
4423914876311586 



In the language of railway transportation, it is equivalent to the fact that the optimal position of junction is 
independent of the currency of the state: it does not matter whether in roubles or in dollars the freight charge is 
nominated. 
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Figure 3: 



1 671 - 79^ / 543856+ 106018^ _ 1 671 - 79^ 671 + 79^ _ 2 
" 3 617 V 83521 ~ 3 617 289 ~ 3 ' 

□ 

Let us discuss now geometrical meaning of the constants from Theorem 15.11 The value ml equals 
the doubled product of the distance from P 1 to P* by the area of the triangle P*P 2 P 3 . 
The first statement of the theorem is equivalent to the fact that 

P *P\ ' ^AP,P 2 P :i + P*P2 ■ S APfP , iPi + PP3 • Sap, Pi p 2 = ^ • 

Finally, the constant ( 15. 2 p is connected with the following one 



h = - 





1 


1 


1 


i 


1 










s 


y* 


y\ 


V2 


2/3 




xl + yl 


xj + yj 


A + vl 


4 + yl 



which is known [13J as the power of the point P* with respect to the circle through the points Pi, P 2 
and P3 (circumscribed circle of the triangle). If one denotes by C the circumcenter of the triangle 
P1P2P3 then 

h = |CP*| 2 - \CPj\ 2 (5.5) 
and provided that P* lies inside this triangle, this value is negative. 
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Results of the present section can evidently be extended to the case of three (and more) dimen- 
sions: 



Theorem 5.2 Let the points {Pj = (xj,yj, Zj)}1j =1 be noncoplanar, and, in addition, be counted in 
such a manner that the value of the determinant 



V 



1111 

X\ %2 X3 Xi 

2/1 2/2 2/3 2/4 

Z\ Z2 z 3 z 4 



(5.6) 



is positive. Then for the choice 



{m* 



\P.Pj\'W j=1 ( 5 - 7 ) 

where Vj equals the determinant obtained on replacing the j-th column of Ii5.6]) by the column^ 
[1, x*, y*, z*] T , the function 

4 

F(x,y,z) = ^m)\J (x - Xj) 2 + (y - y 3 ) 2 + (z - z 3 ) 2 



has its stationary point at 
J5.7}) are all positive and 



3=1 

1 y* •> ^* ) 



If P* lies inside the tetrahedron P1P2P3P4 then the values 







z*) = min F\ 


x,y,z) 




1 


1 


1 


1 


1 






x 2 




X4 


y* 


y-i 


2/2 


2/3 


2/4 




Z\ 


Z2 


Z3 


Z A 


+ vl + 4 


x 2 + y 2 + z\ 


X 2 2/2 ' Z 2 


4 + 2/I + 4 


x^ ~\~ 2/4 z^ 



(5i 



Geometrical meanings of the values appeared in the last theorem are similar to their counterparts 
from Theorem 15 .11 For instance, the value (15.61) equals six times the volume of tetrahedron P1P2P3P4, 
while the value (15.81) divided by V is known [13] as the power of the point P* with respect to a 
sphere circumscribed to that tetrahedron; it is equivalent to (I5.5P where C this time stands for the 
circumcenter of the tetrahedron. 



3 Here T denotes transposition. 
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6 Conclusions 

Analytical solution for the generalized Fermat-Torricelli problem is presented. The three point case is 
completely solved in "extended radicals" : in addition to elementary and extraction of roots operations 
the sign function is utilized in the formulae. For further investigations remains the treatment of the 
multidimensional n point case, although some theoretical results like [lj do not give reason to hope 
for the nice (e.g. extended radicals) solution. 
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